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Abstract 



The paper presents new simple sharp bounds for transition density functions for 
time-homogeneous diffusions processes. The bounds are obtained under mild condi- 
tions on the drift and diffusion coefficients, extending and substantially improving 
PLh ■ previous results in the literature which were limited to drifts satisfying a linear 

growth condition. They lead to an asymptotic expression for the time t transition 
density as t — > 0. While the focus is on the one-dimensional case, an extension to 
multiple dimensions is discussed. Results are illustrated by numerical examples. 
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^". ! 1 Introduction 

o. 

Let iU t ) be a one-dimensional time-homogeneous diffusion process satisfying the stochastic 
differential equation 

XI : dU t = v(U t )dt + <x(U t )dW t , U = u , 

where (W t ) denotes a Brownian motion. The aim of this note is to bound, from above and 
below, the transition probability density function for (U t ), pu(t,Uo,w) := -£^Pu(t,u ,w), 
where Pu(t, Uq, w) := F(U t < w\Uq = Uq). While the focus is on the one-dimensional case, 
the results are easily extended to some special cases in M n , n > 2 (see remark at the end 
of this section). Some simple bounds for the distribution function are also considered. 

Except for a few special cases, the transition functions are unknown for general diffusion 
processes, so finding approximations to them is an important alternative approach. We 
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use Girsanov's theorem and then a transformation of the Radon-Nikodym density of the 
type suggested in pQ to relate probabilities for a general diffusion {U t ) to those of a 
'reference diffusion'. Using a reference diffusion with known transition functions, we are 
able to derive various bounds for the transition functions under mild conditions on the 
original process. The results have a simple form and are readily evaluated. 

As an aside, the generator of the diffusion (U t ) is given by 

*(x) = „(x)§Ui<r'(.)g, 

and the transition probability density function is the minimal fundamental solution to 
the parabolic equation 

A-^)u(t,x)=0. 

Thus the results presented here also bound solutions to certain types of parabolic partial 
differential equations. 

Several papers on this topic are available in the literature, especially for bounding the 
transition probability density. Most recently, [8] proposed upper and lower bounds for 
diffusions whose drift satisfied a linear growth constraint. This appears to be the first 
such paper to relax the assumption of a bounded drift term. The results in [S] will be 
compared with those obtained in the current paper, although the former can not be used 
for processes not satisfying the linear growth constraint. To the best of our knowledge, 
the bounds presented in the current paper are the only ones to relax this constraint, and 
also appear to generally offer a tightening of the bounds previously available. For further 
background on diffusions with bounded drift, see e.g. [H] and references therein. 

In addition, the same ideas allow us to obtain bounds for other functions related to the 
diffusions. This is not the focus of this note and is not discussed in great detail here, but 
as an example at the end of Section [2] we consider the density of the process and its first 
crossing time. This has application in many areas, such as the pricing of financial barrier 
options. Bounds for other probabilities may be derived in the same manner. 

Consider a one- dimensional time-homogeneous non-explosive diffusion (Ut) governed by 
the stochastic differential equation (SDE) 

dU t = v(U t )dt + a(U t )dW u (1) 

where (Wt) is a Brownian motion and cr(y) is differentiable and non-zero inside the dif- 
fusion interval (that is, the the smallest interval I C R such that U t G I a.s.). As is 
well-known, one can transform the process to one with unit diffusion coefficient by letting 

F ( y):= [ y ^s du (2) 

Jyo a \ u ) 



for some yo from the diffusion interval of (Ut) and then considering X t := F(U t ) (see 
e.g. [llj, p. 161). By Ito's formula, (X t ) will have unit diffusion coefficient and a drift 
coefficient fi(y) given by the composition 

^y):=(^--\^oF-\y). 

From here on we work with the transformed diffusion process (X t ) governed by the SDE 

dX t = n(X t )dt + dW u X = F(U ) =: x. 

Conditions mentioned throughout refer to the transformed process (Xt) and its drift 
coefficient fi. 

We will consider the following two cases only (the results extend to diffusions with other 
diffusion intervals with one finite endpoint by employing appropriate transforms): 

[A] The diffusion interval of (Xt) is the whole real line R. 

[B] The diffusion interval of (Xt) is (0, oo). 

The results extend to diffusions with other diffusion intervals with one finite endpoint by 
employing appropriate transforms. 

For the diffusion (Xt) we will need a reference diffusion (Y t ) with certain characteristics. 
The reference diffusion must have the same diffusion interval as (Xt) and a unit diffusion 
coefficient, so that Girsanov's theorem may be applied to (Xt). To be of any practical 
use, the reference process must also have known transition functions. In case [A], we use 
the Brownian motion as the reference process, while in case [B] we use the Bessel process 
of an arbitrary dimension d > 3. 

Recall the definition of the Bessel process (Rt) of dimension d = 3,4,..., starting at 
a point x > 0. This process gives the Euclidean norm of the d- dimensional Brownian 
motion originating at (x, 0, . . . , 0), that is, 



R t = sJ(x + wiy + ...+ (wi d) )\ 

where the (W t ) are independent standard Brownian motions, i = 1, . . . ,d. As is well 
known (see e.g. [10], p. 445), (Rt) satisfies the SDE 

dR t = ^——dt + dW t . 
z ti t 



Note that for non-integer values of d the Bessel process of 'dimension' d is defined using 
the above SDE. The process has the transition density function 

PR (t,y,z) = zfy\-h-^»i%( 1 f 

where 77 = d/2 — 1 and T v (z) is the modified Bessel function of the first kind. For further 
information, see Chapter XI in [10j . 

We denote by W x and K x probabilities and expectations conditional on the process in 
question ((X t ) or some other process, which will be obvious from the context) starting at 
x. We work with the natural filtration T s := a (X u : u < s). 

Finally, note that the present work can be easily extended to a class of n-dimensional 
diffusions for n > 2. If (X t ) is an n-dimensional diffusion satisfying the SDE 

dX t = n(X t )dt + dW t , 

(Wt) being an n-dimensional Brownian motion, then the majority of results can be ex- 
tended assuming fi(-) is curl-free. The extension is straight-forward, and in this note we 
shall only concern ourselves with the one- dimensional case. 

2 Main Results 

This section states and proves a result relating probabilities for the diffusion (X t ) to 
expectations under an appropriate reference measure. In the case [A], the result may 
be known, and we state it here for completeness. The extension to case [B] is straight- 
forward. We then apply this proposition to obtain bounds for transition densities and 
distributions. 

Relation to the Reference Process 

We define the functions G(y) and N(t) as follows, according to the diffusion interval of 
(Xt): 

[A] If the diffusion interval of (Xt) is M, then we define, for some fixed yo G K, 

ry 
G (y) '■= / K z ) dz > 

(3) 

N{t) := [ (i2'(X u )+fi 2 (X u ))du. 
Jo 



[B] If the diffusion interval of (X t ) is (0, oo), then we define, for some fixed d > 3 
(the dimension of the reference Bessel process) and j/o > 0, 

G(y) ■= f (p(z) - ^>j dz, 

N(t) := jT (V(X U ) - {d ~lf 2 ~ 3) + ^Xu)) du. 

Remark 1. For diffusions on (0, oo), the choice of (i is arbitrary subject to d > 3. Therefore 
this choice can be used to optimise any bounds presented in the next subsection. 

Proposition 1. Assume the the drift coefficient jjl of (X t ) is absolutely continuous. Then, 
for any A 6 Tt, 

F X (A) = E x [ e G{X t )-G(x) e -(l/2)N(t) lA ] ^ 

where K x denotes expectation with respect to the law of the reference process. 

Remark 2. In terms of the original process (U t ) defined in ((T]), the condition of absolute 
continuity of fj,(y) requires v{z) and cr'(z) to be absolutely continuous. 

Proof. The proof is a straight-forward application of Girsanov's theorem and its idea is 
similar to the one used in pQ . We present the proof for case [A] , the proof for case [B] is 
completed similarly (see [3] for the general approach). 

Define Q^ to be the reference measure such that under Q x , X = x and 

dX s = dW s , 
for a Q x Brownian motion (W s ). Set 

dw ( r s ~ 1 r s ~) 

Cs '- = W X = 6XP \l ^ X ^ dWu -2j ^(^duj , 

so by Girsanov's theorem under P^ we regain the original process (X s ) satisfying 

dX s = /j,(X s )ds + dW s , 

for a F x Brownian motion (W s ). The regularity conditions allowing this application of 
Girsanov's theorem are satisfied (see e.g. Theorem 7.19 in [7]), since under both F x and 
Q x the process (X s ) is non-explosive and fi(y) is locally bounded so we have, for any 
t> 0, 

P* ( f fi 2 (X s )ds < oo J = Q x ( J n 2 (X s )ds < oo J = 1. 
We then have, under Q x , using Ito's formula and ([3]), 

dG(X s ) = fi(X s )dX s + ^'(X s )(dX s ) 2 

= v(X s )dW s + ^'(X s )ds. (4) 



Note that in order to apply Ito's formula, we only require \x to be absolutely continuous 
with Radon-Nikodym derivative \J (see e.g. Theorem 19.5 in [5]). This also implies the 
above is defined uniquely only up to a set of Lebesgue measure zero, and we are free to 
assign an arbitrary value to /i' at points of discontinuity. 

Rearranging (J3J) gives 

/ fi(X u )dW u = G(X S ) - G(X ) - \ [ fi'(X u )du. 
Jo * Jo 

Hence 

C s = exp |g(X s ) - G(X ) - l -f (fi'(X u ) + fi 2 (X u )) du\ , 

which together with 

completes the proof of the proposition. 

□ 

Bounds for Transition Densities and Distributions 

Define L and M as follows, according to the diffusion interval of (X t ): 

[A] If the diffusion interval of (X t ) is M., then 

L := ess sup (/j,'(y) + /J> 2 (y)) , 
M:=essini(fi'(y)+fi 2 (y)), 

where the essential supremum/infimum is taken over R. 

[B] If the diffusion interval of (Xt) is (0, oo), then, for some fixed d > 3 (the dimen- 
sion of the reference Bessel process), we put 



L := ess sup // (y) — = \- // (y) 



(d- l)(d-3) , u2i 

4y2 



M := ess inf I fi'(y) — 



(d-l)(d-3) 



+ Av) , 



A y2 

where the essential supremum/infimum is taken over (0, oo). 

Note that in what follows, in the case [B], the dimension of the reference Bessel process 
may be chosen so as to optimise the particular bound. Recall also that (Y t ) denotes the 
reference process (the Weiner process in case [A], the <i-dimensional Bessel process in case 

[B]). 



Corollary 1. The transition density of the diffusion (X t ) is bounded according to 

-tL/2 < p X {t,X,w) tM/2 

-e G H- G %(f,x,«;)- ' {) 

Remark 3. The bound is sharp: for a constant drift coefficient /i, equalities hold in (JHJ). 

Proof. Recall (see the proof of Proposition [T]) we only required \i to be absolutely contin- 
uous, and its value on a set of Lebesgue measure zero is irrelevant. Hence L (respectively 
M) gives an upper (lower) bound for the integrand in N(t) for all paths. Applying 
Proposition [I] with A = {X t G [w,w + h)}, h > 0, gives 

inf e Giv) - G( - x) e- tL/2 F x {Y t e [w,w + h)) < F x (X t e [w, w + h)) 

w<y<w+h 

< sup e G{y) - G{x) e- tM/2 F x {Y t e [w, w + h)). 

w<y<w+h 

Taking the limits as h — > gives the required result. D 

In the case of bounded L and M this immediately gives an asymptotic expression for the 
density px(t, x, w) as t — > 0. 

Corollary 2. If — oo < L, M < oo, then, as t — > 0, 

Px (t,x,w)~e G ^- G ^p Y (t,x,w), 
uniformly in x, w. 

While the tightest bounds for the transition distribution are obtained by integrating the 
bounds for the density given above, this does not, in general, yield a simple closed form 
expression. We mention other, less tight bounds that are simple and are obtained by a 
further application of Proposition [TJ 

Corollary 3. The transition distribution function of the diffusion (Xt) admits the follow- 
ing bound: for any toGR, 

inf e G{y) - G{x) e' tL/2 P Y (t,x,w) < P x (t,x,w) < sup e G{y) - G{x) e- tM/2 P Y (t, x, w), 
£<y<w (<y<w 

where £ is the lower bound of the diffusion interval. 

The assertion of the corollary immediately follows from that of Proposition [T] with A = 
{X t <w}. 

By considering other events (e.g. A = {X t > w}), other similar bounds can be derived. 



Further Probabilities 

While the focus of this note is on bounds for the transition functions, Proposition [I] can 
be used to obtain other useful results. For example, consider 

T] X {t, x, y, w) := —F x ( sup X s > y, X t < w ) . 
aw \o< s <t J 

Such a function has applications in many areas, for example the pricing of barrier options 
in financial markets. Using ideas similar to the proof of Corollary [T] immediately gives 

Corollary 4. For the diffusion (X t ), 

e -tL/2 < Vx(t,x,y,w) < e _ tM/2 

- e G( - w ^ G ^r]Y(t,x,y,w) ~ 

Note that for such probabilities the bounds may be improved, if desired, by replacing 
L and M with appropriate constants on a case-by-case basis. For example, if we are 
considering the probability our diffusion stays between two constant boundaries at the 
levels c\ < C2, then the supremum (for L) and infimum (for M) need only be taken over 
the range c\ < y < c 2 . 

Other probabilities may be considered in a similar way. 

3 Numerical Results 

Here we illustrate the precision of the results from the previous section for transition 
densities. Bounds from Corollary [T] are compared with known transition density functions 
and previously available bounds for the Ornstein-Uhlenbeck process in the case [A]. For 
the case [B], we only compare the bounds obtained in the current paper with exact 
results, since there appears to be no other known bounds in the literature. We also 
construct a 'truncated Ornstein-Uhlenbeck' process in order to compare our results with 
other bounds available in the literature. For the Ornstein-Uhlenbeck process we also 
consider an example to illustrate Corollary HI 

The Ornstein-Uhlenbeck Process 

We consider an Ornstein-Uhlenbeck process (St), which satisfies the SDE 

dS t = -S t dt + dW t . 
This process has the transition density 

. . e t I (we 1 — x) 

ps(t,x,w) = exp — — - 

n(e 2t -I) \ 1 - e zt 



see e.g. (1.0.6) in [2J, p. 522, and we begin by comparing this with the bound obtained in 
Corollary [U Since fi(z) = —z, we have 



M 



-1, G(w) - G{x) 



^-w 2 ) 



giving the estimate 



ps(t,x,w) < e2 (x ' 2 w2+t) p w (t,x,w). 



Clearly in this case the bound is tighter for smaller values of \x\ and t. Figure [T] displays 
a plot of the right-hand side of this bound together with the exact density for x = and 
t = 1,2. 





Figure 1: Transition density for an Ornstein-Uhlenbeck process, alongside its upper bound, 
with x = 0. The left-hand side displays the functions for t = 1, the right for t = 2. 



To compare our results with other known bounds for transition functions, we look at the 
bound given by (3.3) in [8] (which, to the best of the author's knowledge, is the only 
bound available for such a process). Figure [2] compares this bound with that obtained in 
Corollary [T] and the exact transition density. The values x = and t = 1 are used (for 
the bound in [8J, q = 1.2 seemed to give the best result, see [S] for further information on 
notation). Note that [8] gives a sharper bound for w close to zero, but quickly grows to 
very large values as |iu| increases, and in general the bounds presented in this note offer 
a large improvement. This is typical for all values of t, with the effect becoming more 
pronounced as t decreases. A meaningful lower bound for this process is unavailable by 
the methods of the present paper, since L = — oo. 

For this example, we briefly look at the bound obtained in Corollary HI We have, see e.g. 
(1.1.8) in [2], p. 522, 



rjs(t,x,0,z) 



V^O" 



-2f 



exp 



(\z\ — xe ') 



t\2 



1-e 



-2t 




0.45 



Figure 2: Transition density for an Ornstein-Uhlenbeck process (solid line) compared to 
bounds given in [8] (dashed line) and Corollary [T] (dotted line), for x = and t = 1. 

Figure [3] compares this as a function of t 6 [0, 1] with the bound obtained in Corollary HI 

1 



T]s(t, x, 0, z) < exp <{ ^(x 2 - z 2 + t) } r] W (t, x, 0, z) 



exp <^ -(x 2 — z 2 + t) 



1 



^vrt 



exp 



1 
2t 



(N-x) 5 



where rjw(t, x, 0, z) is given by (1.1.8) on p. 154 of [2]. 




Figure 3: Bound for r]s(t,x,0,z) compared with its true value, for x = z = —0.5. 

The Truncated Ornstein-Uhlenbeck Process 

Other density bounds available in the literature hold only for processes which have 
bounded drift. For completeness we compare one such bound with the results of this 



10 



paper. We use the bound in [9], which is the most recent for bounded drift and seems 
to give the best results over a large domain. To use these results, however, we need a 
process with bounded drift. As such, we have chosen the 'truncated Ornstein-Uhlenbeck' 
process, which we define as a process (S t ) satisfying the SDE 

dS t = ii(S t )dt + dW t , 

where, for a fixed c > 0, 





c, z < — c, 


fi(z) = 


- < —z, \z\ < c, 




[ — c, z > c. 


For this process we again have M = 


— 1 and, assuming \x\ < c, 




( \[c 2 + x 2 ) + cw, w < — c, 


G(w) - G(x) = < 


\{x 2 — w 2 ), \w\ < c, 




[ \{c 2 + x 2 ) — cw, w > c. 



Figure H] displays the bounds from Corollary [T] together with those in [H] for different 
values of c with x = and t = 1. Smaller values of c move the bounds closer together, 
however for the given choice of x and t they do not touch until we use the (rather severe) 
truncation c ~ 0.45. In general the method outlined in this note provides a dramatic 
improvement. We have also plotted an estimate for the transition density using simulation. 
The simulation was performed using the predictor-corrector method (see e.g. [6] p. 198), 
with 10 5 simulations and 100 time-steps. 

A Diffusion on (0, oo) 

Finally we consider a process from the case [B]. The author believes this is the first paper 
to present a bound on transition densities without the linear growth constraint. The 
process (Vj) satisfying the SDE 



dV t = {pV t + q)dt + ^rV t dW t (6) 

with p, gel and r > 0, has a known transition density (see (26) in [4|). After applying 
the transform Z t = F(V t ), with F(y) = y-y by ([2]), we obtain the process 



dZ t = fjt(Z t )dt + dW t . 
where 

2 y \r 2 

For q > r this dominates the drift of a Bessel process of order 2q/r > 2 so is clearly a 
diffusion on (0, oo). 
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Figure 4: Simulated density and bounds for the transition density of the truncated 

Ornstein-Uhlenbeck process (St)- The solid lines give the simulated densities, the dotted lines 

the bound given in Corollary Q] and the dashed lines the bounds from [9]. Both graphs display 

the functions for x = and t = 1, while the left graph displays them for c = 1, the right for 

c = 2. 



We take the values q = 2.5, r = 1 and p = 1. Using these values, we have 

- 1 (_ (d-l)(d-3Y 



M 



inf 

0<y<oo 



V - + 2.5 + - ( 2 - 
4 y 2 



and 



G(y)-G(x) = \(y 2 -x 2 ) + c\og^- / 

where d is the order of the reference Bessel process and c = 2 — (d — l)/2. 

It remains to choose the order of the reference Bessel process. It is not clear how to define 
the 'best' order of the reference process for a range of z values, as for fixed t and x the 
upper bound for pz(t, x, w) is minimised for different values of d depending on the value 
of z. In Figure E] we have taken t = x = 0.5 and used d = 4.7, however depending on the 
relevant criterion improvements can be made. Again, a meaningful lower bound for this 
process is unavailable by the methods of this paper, since L = — oo. 
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Figure 5: Transition density for the diffusion ([6]), alongside its upper bound, with x = 0.5 

and t = 0.5. 
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